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ABSTRACT 

We use the pulsational properties of the RR Lyrae variables in the globular cluster NGC 1851 to 
obtain detailed constraints of the various sub-stellar populations present along its horizontal branch. 
On the basis of detailed synthetic horizontal branch modeling, we find that minor helium variations 
(Y^O. 248-0. 280) are able to reproduce the observed periods and amplitudes of the RR Lyrae variables, 
as well as the frequency of fundamental and first-overtone RR Lyrae stars. Comparison of number ra- 
tios amongst the blue and red horizontal branch components and the two observed subgiant branches 
also suggest that the RR Lyrae variables originated from the progeny of the bright subgiant branch. 
The RR Lyrae variables with a slightly enhanced helium (Y^O. 270-0. 280) have longer periods at a 
given amplitude, as is seen with Oosterhoff II (OoII) RR Lyrae variables, whereas the RR Lyrae vari- 
ables with F^O. 248-0. 270 have shorter periods, exhibiting properties of Oosterhoff I (Ool) variables. 
This correlation does suggest that the pulsational properties of RR Lyrae stars can be very useful for 
tracing the various sub-populations and can provide suitable constraints on the multiple population 
phenomenon. It appears of great interest to explore whether this conclusion can be generalized to 
other globular clusters hosting multiple populations. 

Subject headings: globular clusters: general — globular clusters: individual(NGC 1851) stars: abun- 
dances, distances, Population II 



1. INTRODUCTION 

Especially following the detection of two distinct sub- 
giant branches (SGBs) in the color-magn itude diagram 
(CMD) of NGC 1851 (TMilone et all 120081 ) . attempts to 
piece together the formation history of this cluster have 
become alluring. One promising explanation of the split 
between the bright SGB (SGBb) and faint SGB (SGBf) 
is that the two subpopulations differ in age by about 1 
Gyr, and t his scenario has been discussed in a numb er 
of studies dMilone et all 120081 iCarretta et "all 120 11 al ibi: 
iGratton et al.ll2012a[ ). Another valid explanation is that 
the SGB splitting is due to differing C+ N+O contents 
and that the two SGBs a re nearly coeval (jCassisi et al.l 
l2008HVentura et al.|[2009h . The horizontal branch (HB) 
of NGC 1851 is also bimodal, with both a prominent 
red HB clump and a blue tail. From the morphol- 
ogy of the HB and the main sequence (MS), strong he- 
lium variations within the cluster do no t seem likely 
(jSalaris et al.l 120081: iD'Antona etall I2009D . and recent 
spectroscopy of the blue HB stars sugge sts minor he- 
lium enhancements (jGratton et al.l l2012af ). Lastly, the 
red giant branch (RGB) is known to harbor different 
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populations (iGrundahl et al.lll999tlCalamida et al.ll2007l : 
iLee et alJlTooa lHan et al.ll2009t ). 

The stellar distribution along the HB of globular clus- 
ters is com monly used to under s tand their formatio n and 
evolution (jGratton et al.l 120101: UTotter et al.ll2010D . and 
previous pape rs dealing with th e modeling of the HB 
of NGC 1851 (jSalaris et al.llifflOl IGratton et al.N2012al ) 
used synthetic HB models to obtain scenarios of the for- 
mation this clusters bimodal horizontal branch. In this 
paper we will approach a more detailed investigation of 
the portion of the HB dealing with the instability strip 
(IS), by discussing the case of the RR Lyrae period dis- 
tribution in NGC 1851. 

Period distributions of RR Lyrae stars have been 
shown to place strong constraints in the framework of 
canonical HB evolution. For example, the problem 
of the peaked distribution of the R R Lyrae periods 
in M 3 (|Castellani fc Tornambel 119811: iRood fc Crocked 
11989( 1 has challeng e d mo del predictions as nicely de- 
scribed by iCatelanl (|2004h . Recent studies have come 
to explain both the M3 period distribution and HB mor- 
phology as the consequence of a range of initial He to- 
gether with a uniform total RGB mass loss (with a very 
small spread) (jCaloi & D'Antonall2008h . or a suitable bi- 
modal mass-loss effici ency along the RGB but a single 
initial He abundance (Cas tellani et al.ll2005| ). 

The RR Lyrae properties of NGC 1851 have been used 
to describe this cluster as "truly an unusual Oosterhoff 
type I object" ([Downes et al.l 120041) . This is largely be- 
cause NGC 1851 is at the extreme end of the Ool-type 
clusters, with its RR Lyrae variables having not only a 
longer than average period for its Oosterhoff class, but 
also a ratio of first overtone RR Lyrae (RR1) to funda- 
mental mode R R Lyrae (RR O) stars more in line with 
OoII- type GCs (|Walkeriri998l) . The suggestion has also 
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been made that the RR Lyrae variables can be divided 
into two subgroups based on their Ca uvby photome- 
try ([Lee et al.l I2009D . although they acknowledge that 
their sample is small, and the apparent bi-modality may 
merely reflect a calcium metallicity spread in the vari- 
ables. Further suggestions that this cluster may be differ- 
ent than other Galactic GCs comes from its phase-space 
distribution, which indic ates it may be associated with 
the C anis Major dwarf (|Frinchabov et al.l [20041: iMartinl 
I2004L however see also Lopez-Corredoira et al. 2007 who 
show that the signatures of the Canis Major dwarf can 
be fully accounted for by Galactic models without new 
substructures). 

Reproducing theoretically both the morphology of its 
unusual CMD as well as the period distribution of its RR 
Lyrae variables is an important step in piecing together 
the formation scenario of NGC 1851. 

2. RR LYRAE OBSERVATIONS 

2.1. Sample and Completeness 

The most complete study of the RR Lyra e variables in 
NGC 1851 was carried out by I Walker! (| 19981) . He presents 
33 variables in a 13.6 arcmin 2 area centered on this clus- 
ter in the B,V and J passbands, 30 of which are RR Lyrae 
stars. Recently ISumerel et al.l (120041) dis covered 19 addi- 
tional variables an dlD ownes et al.l (|2004l) reported eleven 
variables, all within 40" of the cluster center. There is 
overlap in these two samples , as described b y the 2011 
update of NGC 1851 in the iClement et all (j200l cat- 
alog, and most of the new discoveries are RR Lyrae 
stars, although the classification for a handful o f these 
sta rs is still uncertain Neither ISumerel et al.l (2004) 
nor lDownes et all (|2004[) provide calibrated mean magni- 
tudes or amplitudes for their variables, making it difficult 
to use these sta rs in ou r anal ysis. As there is no indi- 
cation that the iWalkerl (|1998| ) sample is incomplete at 
distances greater than 40" from the cluster center, we 
limit our sample of RR Lyrae stars with which to com- 
pare our HB models to this out e r regi on. Of the 29 RR 
Lyrae stars studied by IWalkerl ([1998D that are greater 
than 40" from the center, 25 have both unblended mag- 
nitudes and well determined periods. Thus our sample 
of variables is 25/29 or 86% complete (outside the inner 
core). 

The position on the sky of our sample of variables is 
shown in Figure [TJ and the central 40" is designated by 
a circle . Ligh t curves for these stars are presented by 
IWalkerl (|1998l ) in the BVI passbands, so for these 25 RR 
Lyrae variables, robust mean magnitudes, periods and 
amplitudes are available. The edges of the instability 
strip, judged by the measured colors of variabl es near 
the st rip boundaries, were also determined by IWalkerl 
(|1998|) . as well as the RR1 - RR0 boundary. 

2.2. Period- Amplitude Diagram 

lClement~fc Shclton (1999a) show that the position of 
an RR Lyrae variable in the period-amplitude (PA) di- 
agram is not a function of metal abundance, but rather 
of Oosterhoff type, and derive PA relations for Ool and 
OoII-type RR Lyrae stars empirically. More recently 
iCacciari et al.l (|2005l) study the PA plane of 3 typical 
Ool-type GCs, 3 typical OoII-type GCs and 3 inter- 
mediate types and find that there is a unique period- 
amplitude relation independent of metallicity for RR0 




Fig. 1. — Image of NGC 1851 showing the 25 RR Lyrae variables 
used in our analysis. The central 40" is designated by a circle. 
This region is known to harbor RR Lyrae variables, but due to the 
severe crowding, no calibrated photometry or amplitudes exist for 
the RR Lyrae stars in this region. 

variables in Ool-, OoII- and intermediate-type clusters. 
The periods and V-amplitude of our sample of RR0 
Lyrae variables is shown in Figure [3J and the period- 
amplitude relation of typical Ool and OoII-type systems 
is over-plotted. We note that although many of the RR 
Lyrae stars in NGC 1851 have periods and amplitudes 
that cause them to fall near the Ool PA relation, there 
are a number of stars following the OoII PA relation. 

It is well known that the Blazhko effect, or other ef- 
fects such as a rapidl y changing period, can caus e scat- 
ter in the PA plane ([Clement fc Sheltonl H999al) . The 
Blazhko effect causes the amplitude of light variation to 
vary over t imescales long er than the basic pulsation pe- 
riod. The IWalkerl ([19981) RR Lyrae variables were ob- 
served over an ample time frame (126 total frames ob- 
served over 15 nights during a 1.5 year time span), so de- 
termining the amplitudes using the average light curves 
of the RR Lyrae stars is straightforward. Nevertheless a 
visual determination of the change in amplitudes in each 
RR0 is obtained and shown as an error-bar in Figure [5J 
The change from the average light curve amplitude to 
the maximum Blazhko amplitude ranges from ^0.1 - 0.3 
mag, and photometric uncertainties lead to amplitude 
uncertainties of ~0. 01-0. 05 mag. We conclude that even 
when taking amplitude variations into account, the vari- 
ables are both Ool and OoII-type. 

In comparison, Figure [3] shows the periods and V- 
amplitudes of 1097 RR0 Lyrae variables in 39 Galactic 
globular clusters. The data for this diagram come from 
Table [TJ where each cluster is listed along with the num- 
ber of RR0 Lyrae stars within the cluster that have well 
determined periods an d ^-amplitudes . For complete- 
ness, the [Fe/H] from iCarretta etahl (|2009D . HB-type 
and Oosterhoff type for each GC are also given. This 
sample of RR0 Lyrae variables are divided by their po- 
sition in the period-amplitude plane following the lines 
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Fig. 2. — Top: Period-amplitude diagram for our sample of RRO 
Lyrae variables in NGC 1851. Bottom: Period-amplitude relation 
for 1097 RRO Lyrae variables in 39 Galactic GCs. Our division 
of Ool and OoII-type RR Lyrae variables are shown by dark and 
light ci rcles, respectively. The lines derived by [Clement & Sheltonl 
Ill999al) for Oosterhoff I and Oosterhoff II RRO are overplotted. 

that lClement fc Sheltonl (|1999aD derived for Oosterhoff I 
and Oosterhoff II RRO stars (see Figure [2]). Here an Ool 
RR Lyrae variable is denned by 

A(V) > -5.1453 P + 4.02 (1) 

and an OoII-type RR Lyrae variable by 

A(V) < -5.1453 P + 4.02, (2) 

where A(V) is the ^-amplitude and P is the period. 
We define an Oosterhoff ratio for each GC, which is 
simply the number of Ool-type RRO Lyrae stars com- 
pared to the total number of RRO Lyrae stars in the GC, 
OolRRo/TotRRcv 

We note that iCacciari et all ([20051 ) showed for the M3 
RRO Lyrae stars, a quadratic PA relation is a closer fit 
than a linear one. However, their relation does not ap- 
proximate high amplitude RRO Lyrae variables well (Ay 
> 1.5 mag), largely because such variables are absent in 
M3. As our comprehensive sample includes a handful of 
such stars (one Ool- and three OoII-type RR Lyrae), the 
ICacciari et al.l (|2005l ) relation is not used here. 

Figure [3] shows the histogram of the Oosterhoff ratio 
of the GCs in our sample. Most of these clusters lack a 
complete sample of RR Lyrae variables and at least some 
of the RR Lyrae amplitudes arc likely affected by the 
Blazhko effect or other light curve "noise" . But even with 
these caveats, it is clear from the figure that the Ooster- 
hoff ratio splits the GCs into two groups; the Ool-type 
clusters have RR Lyrae variables with shorter periods 
for a given amplitude and hence have larger Oosterhoff 
ratios (with respect to the OoII-type clusters). Further, 
there is an absence of clusters falling in the "gap". We 
therefore believe that our Oosterhoff ratio is useful to 
distinguish between Ool- and OoII-type GCs. Moreover, 



this ratio may be used to evaluate the degree for which 
a GC is a typical Ool- or OoII-type cluster. 
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Fig. 3. — A histogram of the ratio of Ool stars in the Milky Way 
GCs. The Oosterhoff ratio of the NGC 1851 RRO Lyrae stars is 
highlighted. 

The majority of Ool-type GCs have a RR Lyrae star 
population in which OolRRo/TotRRo > 0.8 (i.e., 80% of 
the variables can be defined by Equation 1). In contrast, 
the Oosterhoff ratio of NGC 1851 is 0.74. Equation 1 
and 2 are defined somewhat arbitrarily, and therefore we 
experimented with a variety cuts (where we modified the 
zero-points and slopes) to distinguish between Ool- and 
OoII-type stars. The Oosterhoff ratio of NGC 1851 var- 
ied between ~0. 63-0. 75, a percentage that is smaller than 
83% of the rest of the Ool-type clusters. This indicates 
that NGC 1851 contains variables with more OoII-like 
periods and 1^-amplitudes than the majority of the other 
Milky Way Ool-type GCs. 

The other three Ool GCs that have comparatively 
small Oosterhoff ratios are NGC 362 (0.71), NGC 4590 
(0.58) and NGC 6362 (0.67). It has been noted that the 
RR Lyrae variables of NGC 362 and NGC 1851 are re- 
markably similar in the period-amplitude diagram, sug- 
gesting similar masses and luminosities. This is also seen 
here, as they have very similar Oosterhoff ratios. Re- 
cently it was found that like NGC 1851, NGC 362 also 
has a split sub-giant branch, although the SGBf compo- 
nent includes only a few percent of the total number of 
SGB stars of the cluster (just ~2-3%, A. Milone private 
communication) . 

New photometry of the NGC 1851 RR Lyrae variables, 
in which periods and amplitudes of the inner RR Lyrae 
stars have been obtained, suggest that NGC 1851 has a 
slightly higher Oosterhoff ratio of ^0.80 (Amigo et al 
2012, in preparation). We therefore suspect that the 
Oosterhoff ratio of NGC 1851, whereas still lower than 
the majority of the Ool-type clusters, is not as abnor- 
mal as the globular clusters NGC 362, NGC 4590 and 
NGC 6362, discussed above. 
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2.3. First Overtone vs. The Fundamental Mode 

The number ratio of first-overtone RR Lyrae stars to 
total RR Lyrae stars, Ni/Ntot, is usually given to quan- 
tify the frequency of the different RR Lyrae pulsators. 
In general OoII-type GCs have about 2-4 times the fre- 
quency of RR1 stars as Ool-type GCs. This is thought to 
be related to the transition temperature between the in- 
stability strip for first overtone pulsation and fundamen- 
tal one. For example, moving the transition from RRO 
to RR1 variables toward lower temperatures (i.e., trans- 
forming fundamental in first overtone pulsators) has the 
twofold effect of increasing the periods of the RRO Lyrae 
population, as well as increasing the relative number of 
first overtones. This is discussed in detail by Castellani, 
Caputo & Castellani (2003). 

NGC 1851 has a large Ni/N to t with respect to most 
Ool-type GCs. This ratio strongly depends on the com- 
pleteness of the sample, and because NGC 1851 is too 
crowded for g round-based observations to resolve, the 
iWalkerl (|1998l f sample i s incomplete a t distances close to 
the core. Assuming the Walker (1998) RR Lyrae sample 
is complete at distances larger than 40" (see Figured]), 
Nx/Ntot = 0.27, one of the largest for Ool-t ype GCs. 
For example, from the [Castella ni et al.l (|2003|) compila- 
tion of 32 clusters with 12 or more pulsators and well 
recognized period and pulsation modes, only 3 of the 17 
Ool GCs have an Ni/N to t greater than 0.27. Recent 
studies have disc overed RR Lyr a e star s closer than 40" 
from the center. iDownes et all (|2004ft find Ni/N^t = 
0.54, and using Clement's catalog, the value of this ratio 
is ^0.40. These ratios are somewhat uncertain, however, 
because of the lack of mode identification for some of 
these newly d i scover ed pulsators. The periods listed in 
iSumerel et al.l (120041) a r e "ten tative" and the periods de- 
rived bv IDownes et al.l (|2004fl do not cover a full period 
for any of their objects. Further, because these stars lack 
amplitudes, their position in a PA diagram also can not 
be used as a diagnostic to identify fundamentals and first 
overtones. Although their is no consensus on the value 
of Ni/Ntot yet, it is clear that this ratio is at least 0.27, 
and likely even larger. A Ni/N to t larger than 0.27 is also 
consistent with the results from Amigo et al. (2012, in 
preparation) . 

3. SYNTHETIC HB MODELING 

Previous s ynthetic HB models fo r NGC 1851 have been 
presented b y Catel an~et al.1 (11998 | ). an d m ore recently by 
ISalaris et al l (120081) .iHan et all (12001) and lGratton et all 
(|2012aD . ISalaris" et al.l (|2008l ) c ompared their simula tions 
to the HST observations by IMilone et ail (|2008l) and 
found two satisfactory scenarios to reproduce the CMD 
of HB stars. In both of these models, the blue HB, 
red HB and variable stars are predicted to come from 
the SGBf, and the stars from the SGBb are confined to 
the red portion of the observed sequence. This infer- 
ence was base d on the number rat io SGBf/SGBb=55:45 
determined in IMilone etal] (|2008l ). ISalaris et afl ()2008l ) 
also found that the initial He abundance of HB stars 
had to be relatively uniform to reproduce the CMD de- 
rived bv IMilone et al.l ((2008). In these data the measured 
magnitudes and colors of the RR Lyrae population are 
at random phases; therefore the portion of the observed 
HB crossing the IS could not be used for detailed con- 



straints on the models. However, they did verify that 
recent theoretical pulsation models of RR Lyrae stars 
(|Di Criscienzo et all l2004al) predict an instability strip 
for NGC 1851 too red by ~0.03-0.04 mag in (F606W- 
F 8UW) compared to the IMilone et al.l (j2008f ) data. 

IHan et all ((2009), on the other hand, found from their 
UVI photometry that the RR Lyrae variables could 
come from both the SGBb and the SGBf. They con- 
struct two population models for NGC 1851. In the first, 
the second generation population is more enhanced in 
metallicity but not in helium (AZ=0.0004 and Aage = 
O.lGyr), and in the second, both metal and helium abun- 
dances are enhanced (AZ=0.0004, AY=0.05, and Aage 
= O.lGyr). They find that their AZ-only model is in 
conflict with the observed CMDs of NGC 1851, but that 
their AZ+AY model is in good agreement with the ob- 
servations from the MS to the HB. In this model, the RR 
Lyrae variables exhibit different He abundances. Because 
their RR Lyrae variable sample was found at a random 
phase of pulsation, the RR Lyrae colors could not be 
used as a stringent constraint in their models and the 
RR Lyrae star p eriods were not discuss ed. 

More recently iGratton et al.l (|2012al ) c onsidered the 
revise d SGBb/SGBf ratio determined by IMilone et all 
(2009) and introduced new spectroscopic constraints; 
they find a small difference in the iron content between 
the SGBb and SGBf, and argue that an age spread of 
~1.5 Gyr is the most viable explanation for the splitting 
SGB. They also find that the RHB stars separate into 
two groups depending on their O and Na abundances, 
and that the BHB stars are slightly helium enriched as 
compared to the RHB stars. Hence to satisfy these con- 
straints, the HB is modeled with four different compo- 
nents, with the IS originating from the SGBb. That each 
SGB hosts multiple generations of stars is shown from 
spect roscopy of stars on the double SGB (|Lardo et all 
I2012D , making it likely that multiple components may be 
needed to model the HB. 

Our own synthetic HB calculations described below are 
aimed at answering the following question: what is the 
most straight-forward way to reproduce the RR Lyrae 
instability strip of NGC 1851 - and in particular the pul- 
sational properties of its RR Lyrae variables? We seek to 
provide a simple and attractive explanation for the clus- 
ter HB and IS morphology, keeping the number of free 
parameters to a minimum, yet still reproducing the RR 
Lyrae star properties that make this cluster stand out as 
having an unusual Oostcrhoff type. 

3.1. Synthetic HB models 

The HB evolutionary tracks used here are fr o m the 
BaST I stellar library (jPietrinferni et al.l I2004L I2006L 
2009) and have already been d escribed by ISalaris et all 
(2008) and ICassisi et al.l (120081 ). They are also the same 
that lGratton et al.l (|2012a|) employed. Briefly, evolution- 
ary tracks are for a normal a-enhanced ([a/Fe]=0.4) 
metal mixture, with [Fe/H]=-1.31 dex and Y=0.248. 
The HB tracks were interpolated among the models with 
Y=0.248 and additional BaSTI models with Y=0.300, 
to determine HB tracks for intermediate values of Y, 
at the same iron content. Similarly, to determine HB 
tracks with a milder CNO-enhancement, an interpola- 
tion between the reference set and the models with the 
CNO sum enhanced by 0.3 dex ( Pictrinfe rni et aLll2009D 
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is used for a portion of the synthesis, as in lGratton et al.l 
(|2012af) . Hence, the [Fe/H], a-enhancement and CNO- 
cnhan cement is consistent wit h spectroscopic result s 
from lCarretta et al.l (|2011al ) and IGratton et al.l (|2012ah . 
We wish to remind the reader that, as long as the CNO 
sum is unchanged, the effect of the observed CNONa 
anticorrelations (overimposed to a standard a-enhanced 
metal mixture) on the evolutionary tracks and isochrones 
is negligible, and standard a-enhanced models are ade- 
quate to represent the whole cluster population. Only an 
enhancement of the CNO sum requires the calculation of 
appropriate models. On the other hand CNONa anti- 
correlations even at constant CN O affect the bolometric 
corrections of filters like B and U (jSbordone et al.ll2011l) 
at low effective temperatures, but not longer wavelength 

filters. 

The four HB components described by IGratton et al.l 
(2012a) are used as a starting points for our calcula- 
tions. Objects from our synthetic HB that fall within 
the observed IS from iWalkerl (|1998) are considered RR 
Lyrae variables (this region is labeled in Figure and 
their peri od is calculated from the p ulsation equation 
given by IDi Criscienzo et "all (|2004aT). The inten sity 



mean magnitudes and colors given bv I Walked (fl99ll are 
used as a comparison to the synthetic HB, because the 
static magnitudes and colors from stellar evolution mod- 
els a re represented better by in tensity-averaged quanti- 
ties (|Di Criscienzo et al.l I2004al) . Although the IWalkeH 
( 1998) observations include the BVI passbands, we em- 
ploy only the V and I magnitudes, because - as discussed 
before - they are not affected by the observed CNONa 
abundance anticorrelations. 

In addition to the observed V and (V — I) distribution 
of the HB stars in the CMD, and the observed (B :V :R) 
(blue variable :red HB) ratio of (B : V :R) = (33±8:10± 
5:56±11) (in line w ith the results bv lCatelan et alj|1998t 
iSaviane et al.lll998f ). we impose as a further constraint on 
our simulation the observed distribution of the RR Lyrae 
periods. 

As in IGratton et aLl (l2012aD we adopted E(B - 
V)=0.02 mag ( Walkerlll998|) and fixed the apparent dis- 
tance modulus to (m — M)v=15.56 mag by matching the 
observed mean magnitude of the RHB with our synthetic 
counterpart. The (B :V:R) ratio of our 'best fit' simula- 
tion is (27 :9 :64), consistent with the observed value. For 
reasons that will become clear in the discussion that fol- 
lows, we consider a preliminary reference age of 11 Gyr 
for the progenitors of the RHB stars. This implies, for 
the assumed metallicity and a 'normal' y=0.248, an ini- 
tial mass of O.86M0 for the stars at the tip of the RGB. 
The HB co mponents are d e scribed below. 

(1) As in IGratton et al.l ()2012aD . the majority of the 
RHB population is modeled with normal CNO abun- 
dance, a normal y=0.248, and a Gaussian mass distribu- 
tion with < M > = 0.67 ± 0.005 M©. This corresponds 
to a total mean mass loss AM=0.19 M© along the RGB, 
for the assumed 11 Gyr age. 

(2) A smaller RHB subpopulation, that is Ba-rich, 
mak es up —10% of the HB population. It is modelled, 
as in lGratton et al.l (|2012al ). employing a Gaussian mass 
distribution with < M > = 0.65 ± 0.004 M (corre- 
sponding to AM -0.21 Mq) an 0.15 dex enhanced CNO 
abundance, and normal F=0.248. If we assume that the 
mean total mass loss has to be constant among all cluster 



RGB stars - and equal to AM=0.19 M© as determined 
for the rest of the RHB component - the mean value 
of the mass for this HB sub-population implies an age 
— 1 Gyr older for the progenitors of this HB component. 

(3) The horizontal part of the BHB, including the RR 
Lyrae instability strip, makes up —10% of the cluster stel- 
lar content. This component i s the focus here, and is the 
only one that is modified from IGratton et all (|2012aT ). In 
particular, instead of adopting a constant He abundance 
F=0.265, the helium content for stars between the blue 
end of the red clump and the beginning of the BHB blue 
tail has a continuous distribution between F=0.248 and 
0.280. A simple, flat probability distribution for Y and 
constant AM=0.19 M Q (for an age of 11 Gyr) with a 
lcr Gaussian spread of 0.005 M© - as for the RHB stars 
- for all RGB progenitors provide a good match to the 
observed RR Lyrae periods, as discussed below. 

The mean He abundance in the IS is < Y > =0.271, 
clo se to the con s tant a bundance Y =0.265 employed 
by IGratton et al.l (|2012aD for this component, and the 
mean mass is < M > = 0.634 Mq. It is worth noting 
that the observed HB distribution of these stars is well 
matched by both this simulation and the simulation by 
IGratton et al.l (|2012al ) . This spread in He content is nec- 
essary to reproduce the observed period distribution. 

(4) The blue tail of the BHB population makes up 
-20% of the HB stars. As in IGratton et "all (|2012ab 
this component is modeled with normal CNO abun- 
dance and y=0.28. The Gaussian mass distribution has 

< M >=0.59±0.005 M , that for an age of 11 Gyr would 
correspond to a mean AM=0.22 Mq. If the mean total 
mass loss is instead fixed at AM=0.19 M©, this value of 

< M > implies that the progenitors of the BHB blue tail 
stars are —1.5 Gyr older that the RHB component with 
normal CNO and Y. Notice that the constraint on the 
progenitor Y is weak er for BHB stars (see discussion in 
IGratton et al.ll2012al ) and a small spread of order 0.01 
may be present. 

These results from the HB synthetic modeling can be 
interpreted in terms of the progeny of the SGBb and 
SGBf subpopulations (the ratios — 2/3 and — 1/3 of the 
total SGB population, respectively, are adopted as de- 
termined by Milone et al. 2009) as follows: 

(1) The sum of the fraction of stars in the blue tail of 
the BHB and in the mildly CNO-enhanced Ba-rich RHB 
component is —35% of the total HB population. If we 
consider as a reference the CMD of 11 Gyr old SGB stars 
with 'normal' Y and CNO abundance, the progenitors of 
these two HB components will be distributed along a 
fainter SGB than the reference one. In the case of the 
progenitors of the BHB component this is an age effect, 
for a change of Y does not have a major effect on the 
SGB luminosity. For the Ba-rich RHB progenitors the 
reason is the slightly higher age and the mildly enhanced 
CNO abundance, that act both in the direction of pro- 
ducing a fainter SGB. As a result, both the Ba-rich RHB 
stars and the blue tail HB progenitors display an approx- 
imately coincident SGB, that we tentatively identify as 
the SGBf in the cluster CMD. 

(2) The sum of the fraction of stars in the horizontal 
BHB (including the IS) and the RHB stars with normal 
composition amounts to —65% of the total HB popula- 
tion. We identify their progenitors as the stars harbored 
by the SGBb in the cluster CMD. 
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To conclude this section, it is worth noting that the 
exact value of the assumed reference age (11 Gyr) is not 
critical. Had a different age been assumed, i.e. 10 or 
12 Gyr, the previous conclusions will still be valid. The 
only difference is that all values of AM would need to be 
shifted downwards (or upwards) by ^0.02 M Q - to keep 
the mass distribution along the HB unchanged - but the 
interpretation of the results would be identical. Finally, 
the different chemical composition - and small age dif- 
ferences - assigned to the SGBb and SGBf populations 
do not affect substantially SGB and RGB timescales; as 
a consequence, the number ratio SGBb/SGBf will be 
approximately equal to the number ratio of their HB 
progeny. 

3.2. Comparison With Observations 

This synthetic HB model is shown compared to the 
observed one in Figure 0] where an observational scat- 
ter of ay, i = 0.01 mag is assumed. The four com- 
ponents are high-lighted for clarity and the RR Lyrae 
region is labeled. Our full synthetic HB model repro- 
duces t wo peculiaritie s in th e GMD of NGC 1851 pointed 
out by IBrocato et al.l (|1999l ) , namely the clump of stars 
near the red edge of the HB and the slightly tilted HB 
(Ay^O.l mag). Features such as these are present also 
in NGC 6362 (which has an RR Lyrae population with 
properties very similar to NGC 1851), and in the extreme 
cases of NGC 6388 and NGC 6441 (A^*~0.5 mag). 

The focus here concerns the component that includes 
the instability strip, which comes from the SGBb. As 
iWalkerl (|1998t ) mentioned, the ZAHB is very cleanly de- 
fined and is not horizontal, being slightly brighter at 
bluer colors. This is reproduced in our synthetic HB 
model by stars that range in helium abundance from 
F=0.25-0.28, and range in mass from 0.61 to 0.65M Q . 

The theoretical periods and pulsation amplitudes from 
the RR Lyrae variables in our synthetic HB model are 
compared to the observed periods and amplitudes in Fig- 
ure [5] Here the observed periods co me from 2 8 RRO and 
18 RR1 variables as determined by [Walkerl (119981) and 
ISumerel et al.l (|2004[) to encompass all the data avail- 
a ble in the literature ( see the 2011 update of NGC 1851 
in lClement et al.ll2001l) . We also compare the theoretical 
periods with 27 RRO and 18 RR1 variables as determined 
by Amigo et al. (2012, in preparation). These authors 
derive periods for the recently identified inner RR Lyrae 
variables based on light curves with ^200-300 points 
in each of the B- V- and I-passbands. The o b served 
amplitudes come excl usively from the IWalkerl ([1998D 
RRO Lyrae sample, as ISumerel et al.l (|2004l ) present in- 
strumental magnitudes only. The amplitudes from the 
Amigo et al. (2012, in preparation) sample are not 
used, because unlike when determining periods, ampli- 
tudes can_j2e_ajfected_by crowding and blending issues 
(e.g., iMaiaess et al.|[2012ti and we do not have a feel for 
how/if blending affects their (preliminary) amplitude de- 
t erminations. 

IMarconi et al.l (|2011D provided a detailed comparison 
between the impact of the He abundance on the pulsa- 
tion properties of RR Lyrae stars and concluded that He 
content marginally affects the pulsation behavior of RR 
Lyrae stars. They noted that the increase in the aver- 
age pulsational period associated with the He increase 
is only due to the brighter luminosities which character- 



ize He-enhanced evolutionary models. As a consequence, 
from a theoretical point of view the impact of an He- 
enhancement on the pulsation properties of RR Lyrae 
stars can be directly taken into account when adopting 
evolutionary tracks for the appropriate He abundance 
and pulsational model predictions obtained for a canon- 
ical He abundance. 

The periods of the RR1 variables are fundamentalized 
via log Po ~ logP -I- 0.127, where Po is the fundamental 
mode period, and the theoretical periods are calculated 
for all HB stars fall i ng with in the observed IS using the 
IDi Criscienzo et "ahl (|200 4a') RR Lyrae pulsation models. 
The mean observ ed RRO Lyrae p eriod in NGC 1851 is 

< P > = 0.571 (jCatelanll200"9al ). and fro m our model 

< Po > = 0.569 ± 0.006 d is calculated, where the un- 
certainty is the error in the mean. The number of RR 
Lyrae stars in the simulations is ~5 times larger than the 
number of observed RR Lyrae stars. In this way, in the 
synthetic HB model, the effect of statistical fluctuations 
in the number of objects at a given magnitude and color 
is minimized. 

The observed minimum fundamentalized period of the 
first overtones and the minimum fundamentalized period 
from the synthetic HB is very similar, Pp^0.34 days. It 
has been shown that the minimum fundamentalized pe- 
riod is a robust observable to constrain the evo lutionary 
properties of RR Lyrae stars (|Bono "etaIIH99l . so it is 
especially encouraging that the observed and theoretical 
values agree. 

We have performed a Kolmogorov-Smirnov (KS) test 
to establish whether one can reject the null hypothesis 
that the two samples of periods (observed and synthetic) 
come from the same distribution. From a comparison 
betwe en the combined IWalkerl (|1998l) and ISumerel et al.l 
(2004) sample and the synthetic RR Lyrae stars, the KS- 
test returns a probability P=0.86, well above the default 
threshold Pth=0.05 below which one rejects the null hy- 
pothesis. When the Amigo et al. period distribution 
is compared with our theoretical one, if we assume that 
both samples are drawn from the same parent popula- 
tion, the samples will differ by this amount or more 94% 
of the time. On this basis, we find that the synthetic 
periods from our simulated HB and the observed periods 
agree well with each other. 

Theor etical pulsation amp l itudes are also determined 
from the IDi Cr iscienzo et al.l ()2004aD relations, although 
three points should be taken into account: 

1. Theoretical amplitudes depend on the the mixing- 
length, I /H p , (where l/H p is the pressure scale-height), 
which is uncertai n and may change fr om the blue to the 
red edge of the IS (jMarconi et al.l2003| ). Pulsation ampli- 
tudes are affected by I /H p in the sense that a larger value 
of I /Hp translates to smaller pulsation amplitudes as a 
consequence of the larger efficiency of convective motions 
and, hence larger quenching to the pulsation mechanism 
provided by convection. 

2. The synthetic pulsational amplitude - per iod pre- 
dictions provided by IDi Criscienzo et al.l (|2004aD show a 
large scatter, of the order of ^0.2 mag (standard de- 
viation) in the T^-amplitude at a fixed period (see Fig. 
6 of their paper, where the dashed lines represent the 
standard deviation). Since our predictions of the V- 
amplitude in our synthetic HBs are based on an A(V) 
- period relation obtained by an interpolation on the 



Oosterhoff effect in NGC 1851 



7 



15.5 - 



16 - 



16.5 
> 



17 - 



17.5 - 



18 - 



n 1 — r 



T 



~\ — i r 



T 



n 1 — r 



~\ — i r 



n 1 r 



observed 



■ ^ • ^ • 



• 




IS 



° X * X 

X 

°n x ** 




SYNTHETIC 

SGBf 

o 0.590±0.005, 



Y=0.280 



SGBb 

x 0.630±0.006, Y 
* 0.670±0.005, Y 
0.650±0.004, Y 



1 



1 



0.248-0.280 
0.248 
0.248, 

[CN0/Fe] = 0.15 
j i I i i i I 



-0.2 







0.2 



0.4 

V-I 



0.6 



0.8 



Fig. 4. — Qualitative comparison of the observed HB and the syntheti c HB, assuming four separate HB components. The R R Lyrae 
instability strip is marked, and corresponds to that found by the Walker (1998) study of RR Lyrae variables in NGC 1851. The [Walker 
(1998) RR Lyrae variables are designated by large circles. 



data provided bv lDi Criscienzo et al.l (|2004aD . we do not 
expect a great accuracy in our A( V) estimates. In ad- 
dition, iDi Criscienzo et"aT1 (|2004a[ ) have shown that for 
P < 0.68 d, pulsation model predictions for the V- 
amplitude are still more affected by a change in the 
adopted mixing length value (see lower panel of the 
quoted figure): for increasing mixing length values they 
predict a significant decrease in the dependence of the 
^-amplitude on the pulsation period. 

3. When dealing with pulsating structures, the static 
magnitudes differ from the measured magnitudes, which 
are mean quantities averaged over the pulsation cycle. 
In finding theoretical amplitudes, a correction between 
static and intensity-averaged magnitudes is required. 
The discrepancy between static and mean values is a 
function of the pulsation amplitude, and th e corrections 
adopted here come from iBono et al.l (|1995| ). 



When using a l/H p =2.0, the Walker (1998]) V- 
amplitude distribution is similar to that observed, al- 
though the theoretical amplitudes appear to be ~0.1 
mag larger. This is not completely surprising con- 
sidering the scatter in the IDi Criscienzo et "all (|2004a[ ) 
^-amplitude relation as well as the magnitude correc- 
tions discussed above. Extrapolating linearly between 
l/Hp—1.5 and 2.0, an increase in the mixing length 
of 0.1 would cause a decrease of the theoretical am- 
plitudes by ^0.08 mag. Such a decrease would pro- 
vide a satisfactory agreement between theory and ob- 
servations. A 0.1 change in l/H p is well within the 
uncertainties in the mixing length calibration, and not 
nearly l arge enough to affect the predicted pulsation 
period s (IBono &; Stellingwerj Il994t IDi Criscienzo et al.l 
120041* iMarconi fc DeglTnnocentil 120071) . An I /H„=2.0 
was also used bv IBono et al.l (|2007| ) to derive a visual 
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Fig. 5. — A comparison between the theoretical periods and am- 
plitudes from our synthetic HB versus the observed periods and 
amplitudes. A KS test indicates that statistical significance is de- 
tected between the observed and theoretical periods and ampli- 
tudes. 

distance modulus from nonlinear convective models of 
RR Lyrae stars to NGC 1851 of [m - M)y=15.52±0.11 
mag, which is similar to the distance modu lus adopted 
here. In contrast, when using a l/H p =1.5, iBono et al.l 
(|2007h find a distance modulus of (m- M )y=15.40±0.12. 
Therefore we conclude that using a larger value of l/H p , 
(1/ _ffp~2.0), provides a consistent comparison between 
our synthetic HB and pulsational predictions. 

In general, the RR Lyrae variables with Y < 0.27 fall 
in the Ool area of the PA diagram, whereas the RR Lyrae 
variables with Y > 0.27 fall close to the OoII line. As- 
suming that the period-amplitude diagram can be effec- 
tively used to classify RR Lyrae stars into an Ooster- 
hoff type, this means that He and Oosterhoff type are 
correlated in this cluster. This is not completely unex- 
pected, as an increase in He makes RR Lyrae variables 
brighter and, as a consequence, higher helium abundance 
makes their p ulsational period longer (|Bono 
lMarconill2009l ). 

On a general ground, for a given total mass the HB 
stars with Y"<0.27 are redder than those with Y>0.27. 
The red part of the IS, where the fundamental mode RR 
Lyrae stars reside, is consequently more populated, and 
a smaller Ni/N to t is obtained. Our synthetic HB yields 
Ai/At t ~ 0.1 for the stars with Y <0.27. This is a ratio 
that is seen for the majority of the Ool-type GCs. In 
contrast, the blue part of the horizontal branch contains 
more stars with higher helium abundances, and as the 
first overtone RR Lyrae variables reside in the blue part 
of the IS, a larger Ni/N tot is obtained. Our synthetic 
HB yields AV^tot^O^ for stars with Y>0.27, a ratio 
more in line with OoII-type GCs. The observed ratio of 
first overtone to total RR Lyrae variables for NGC 1851 
in our sample is Ni/N tot ~0.30 and is easily explained 
(and reproduced with our synthetic HB) by the spread 
in Y along the RR Lyrae instability strip. 



We find that simulations using a constant helium 
for the portion of the HB containing the IS (as in 
iGratton et al.ll2012af ) do not fit the constraints given by 
the NGC 1851 RR Lyrae variables as well. For example, 
adopting Y=0.265 results in an Ni/N^ = 0.11 (versus 
the observed Ni/ N to t = 0.28) 

We note that iMilone et all (|200&) provide an up- 
per limit to a possible dispersion in helium abundance 
of AY=0.026 between the two SGBs in NGC 1851, 
a value close to the spread assumed in the syn- 
thetic HB presented here. Other estimates of the He 
spread in NGC1851 are slightly l arger, i.e., AY= 0.04 
(I Ventura et al.1 [20091). AY =0.05 (jHan et all [20091 ) or 
AY=0.048 (|Gratton et alllMoh . 



4. DISCUSSION AND CONCLUSIONS 

The population distribution of the stars along the HB 
has been modeled assuming the presence of (at least) 
four populations with differing helium contents. In our 
simulations, the only parameter we vary is the initial He 
abundance of the HB progenitors, keeping the same total 
RGB mass loss for all components. Both the RR Lyrae 
period distribution as well as the number ratio of first 
overtone RR Lyrae to total RR Lyrae stars, Ni/N to t, 
provides constraints pertaining to the component of the 
HB containing the IS. It is straight forward to repro- 
duce the observed distribution of RR Lyrae stars inside 
the instability strip with minor He variations (Y^O.248- 
0.280) and from a HB subpopulation corresponding to 
the progeny of a fraction of the SGBb stellar population. 

Therefore, the IS of NGC 1851 belongs to a second 
generation (SG) of stars. That a SG exists within 
the SGBb component is in agreement with recent spec- 
troscopy of the SGB stars, showing that each SGB host s 
multiple subgener a tions of stars (jLardo et al.1 I2012H . 
iD'Antona fc Calol (|2008l ) have also postulated that the 
longer periods of the NGC 1851 RR Lyrae may indicate 
that these variables may belong to the SG, and our syn- 
thetic HB strengthens this notion. Chemical anomalies 
in GCs suggest that self-enrichment is a common fea- 
ture among GCs. The quasi-constancy of heavy metals 
in most GCs leads to the assumption that abundance 
variations are not or scarcely affected by SN ejecta, but 
involves formation of a second generation of stars from 
matter processed into the FG stars. 

The SG will most pr obably show a spread in He 
(jD'Antona fc C aloi 2001) because the self-enriched ma- 
terial may come from different progenitors that have dif- 
ferent chemical peculiarities, or may be diluted in differ- 
ent fractions with matter from the first generation (FG) . 
We remark that such a helium spread is an essential in- 
gredient in order to reproduce the pulsation properties of 
the RR Lyrae population as a whole. Simulations using 
constant He across the IS give synthetic period distribu- 
tions that do not match the observed one as well and 
result in the observed Ni/N to t being lower than what is 
observed. Actually, one can note that variations of He 
in 19 GCs have also been deduced by iBragaali a et al.l 
(pOlOh fr om 1400 RGB stars with abundance determi- 
nations. As discussed in IGratton et al.1 (|2010l ). a star- 
to-star spread in the He abundance may explain many 
aspects of the horizontal branches of GCs. 

It is worth pointing out that there have been sugges- 
tions of problems in the late stages of HB evolution in 
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curre nt HB trac ks dCat elan 2009bT; iValcarce fc Catelanl 
120081) . Moreover, [Catelasl (l2009^show that in the case 
of NGC 1851, over- luminous stars on the blue HB could 
be interpreted by an underestimate of the luminosity evo- 
lution along the HB rather than in terms of a moder- 
ate level of helium enrichment. Here we do not attempt 
to resolve this ambiguity for NGC 1851; rather, we as- 
sume that the evolutionary tracks adopted represent the 
HB evolution accurately, and remind the reader that our 
comparisons are ultimately subject to both theoretical 
and observational uncertainties. 

The pulsation periods and amplitudes from the RR 
Lyrae variables resulting from variations in He along the 
IS have different characteristics. The RR Lyrae variables 
with a "normal" helium have periods and amplitudes, 
as well as a Ni/N to t ratio, that is inline with Ool-type 
GCs. In contrast, the RR Lyrae variables with slightly 
enhanced He (0.27 < Y < 0.28) have longer periods and 
a higher ratio of Ni/N tot , indicative of RR Lyrae vari- 
ables in OoII-type GCs. In the absence of spectroscopy 
of the RR Lyrae variables in NGC 1851, the synthetic 
horizontal part of the HB and RR Lyrae instability strip 
presented here is the simplest one that reproduces the 
available observations with the smallest amount of free 
parameters. New observations of the RR Lyrae variables 
may require more complex modeling, however, and would 
be particularly interesting. 

Oosterhoff-I clusters tend to be more metal-rich and 
host fainter RR L yrae variables than OoII clusters 
(|Caputo et al.|[2000l ). As the metallicity has an effect 
on the absolute magnitude of an RR Lyrae, it has been 
difficult to disentangle whether the metallicity difference 
alone is affecting the brightness differences, or whether 
there are differences in the intrinsic magnitudes of RR 
Lyrae variables in Ool and OoII globular clusters caused 
by something other than just metallicity (like evolution 
or helium). In this cluster, where an internal spread in 
[Fe/H] is small at most, our results indicate that a dif- 
ference in helium abundance in the RR Lyrae variables 
is affecting where in the PA diagram the RR Lyrae star 
falls. 

Our model consists of a BHB that is He-enriched 
(Y -0.28) yet older than the RHB (F=0.248). This 
can be explained if the cluster was formed by a process 
such as a merger with populations that differ in He and 
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age. Such a scenario h as a lready been discussed by e.g. 
iCarretta et all (|2011bD and I Bekki feTonel (|2012T ). Hence 
the BHB would not be a second generation (SG) of stars 
originating from the same population as the RHB (the 
SGBb). 

Since this paper was submitted, results from an in- 
termedia te resolution spectrosc opic analysis of the two 
SGBs bv lGratton et all (|2012bH indicate that the [Fe/H] 
difference between the SGBb and SGBf is —0.07 dex, the 
SGBf being more metal rich. The RR Lyrae stars in our 
scenario are the progeny of SGBb; hence this metallicity 
difference does not affect our results. We find that for 
BHB stars hotter than the IS a 0.1 dex increase in [Fe/H] 
at fixed (V — I) changes the HB masses (at fixed Y) by 
— 0.01M Q , and My changes by —0.01 mag. Therefore the 
effect of such a [Fc/H] difference between the two SGBs 
has a negligible effect on our HB modeling. 

We have shown that a spread in He reproduces the pul- 
sational properties of the RR Lyrae sample as a whole, 
indicating the presence of a SG of stars in NGC 1851. 
Our analysis therefore demonstrates that RR Lyrae prop- 
erties in a given GC can provide suitable constraints on 
the multiple population phenomenon in that GC. It is 
worth carrying out more studies of this kind to investi- 
gate further this connection with the occurrence of the 
multiple population phenomenon, especially in GCs with 
a sizable population of RR Lyrae stars and in which the 
stellar chemical patterns are well known. 
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Oosterhoff effect in NGC 1851 



TABLE 1 

MW GCS WITH A POPULATION OF RR LYRAE VARIABLES 
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Oo-Typo 


N R] 
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-0 


.87 
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M68 


-2 


.27 


+0 


.17 


OoII 


12 


IC 4499 


_ 


-1 


.62 


+0 


11 
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62 


Ruprccht 106 
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-1 


.78 
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.30 
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